5922 


Copyright © 1995 by the 




Theoretical Study of Fe(CO)„ 


NASA-TM-1118A0 


Alessandra Ricca and Charles W. Bauschlicher, Jr.* 

NASA Ames Research Center, Moffett Field. California 94035 
Received: October 19, 1994: In Final Form: February 2, 1995 ® 



The structures and CO binding energies Fe to CO 

(DFT) approach. The structures and groun s trends in the CO binding energies for Fe(CO)„ 

2** donation and minimizing Fe-CO 5 a repuls.on. The *endsm £ COb^ ^ ^ 

and the differences between the trends for Fe(CO „ '^enl with experiment, while the first is too 

a perturbational estimate of the connected triple excitations. 

Leopold 9 have obtained the negative ion photoelectron spectra 
of FeCO - at an instrumental resolution considerably improve 
ier that employed by Engelking and Lineberger. They have 
assigned the ground state of FeCO as a I state wi an 

occupation of 3 dd Wsda 2 spa' and meyhave^gge^ted *at 

the bonding is derived from the 4 F state of Fe . They measured 
the Fe-C stretching frequency (465 ± 10 cnT'land the _bmbng 
frequency (230 ± 40 cm-) of FeCO'. ^ jujh ^C 
stretching frequency is consistent with a strong (33.7 ± 3.5 kcal/ 
mol 3 ) Fe-CO bond. 

Castro Salahub, and Fournier 1 ‘have studied FeCO using a 

linear combination of Gaussian-type orbital-density functional 
(LCGTO-DF) method. Their results confirm the assignment 
by Villalta and Leopold of a ground state. Castro « al. 
obtain an Fe-C stretching harmonic frequency (566 cm ) 
higher than the experimental fundamental and a bending 

frequency (272 cm” 1 ) close to that from experiment. Their CO 
binding energy was in good agreement with experiment 
However, it should be noted that when they applied ihesame 
method to FeCO + , they obtained a binding energy th 
significantly too large. 


I. Introduction 

The successive ligand binding energies offer insight into how 
the bonding in metal-ligand systems is changing with the 
number of ligands. For example, a large change in the binding 
energies could indicate a change in the bonding mechamOTV 
For Cu(H 20)„ + there is a large decrease for the ir wa 
binding energy. Calculations have shown that 
tion which is very important in reducing the metal ligand 
repulsion for the first two waters, is lost with the addition o 
the third ligand and is therefore responsible for the large decrease 
binding nnnrgg. For «CO>,* Jta. .s a snjn to 
drop 2 in binding energy for the third CO. Because Fe( )z 
is a quartet state and Fe(CO)s + is a doublet state it isdifficult 
to identify the origin of the decrease m the third ligand binding 
energy for Fe(CO)„ + ; it could be from the loss of so 
Sization as in Cu(H 2 0)„ + or it could be due to a change n 
the Fe spin state. Calculations have shown- that the decrease 
hTthe third CO binding energy in Fe(CO)„+ is due to loss of 
sda hybridization and that the decrease m the fifth CO binding 
enercv is due to the change in the Fe spin. 

These examples show the importance of performing accurate 
calculations to aid in the understanding of the changes in 
bonding that occur as the number of ligands is changed- 
Conversely it is important to have accurate experimental binding 
energies to compare with the calculations. Agreement between 
theory and experiment shows that theory is correctly descnbmg 

^stcessive CO bond energies of Fe(CO)r have been 
determined experimentally by Sunderlin, Wang ajd Squires^ 
The CO binding energies increase forthefrst toee 
molecules and then decrease slightly for the fourth CO. Thi 
^d i very different from that observed for Positive ions 

“libg F«. Stob * F«CO),V-‘ 

(CHA, + 8 where there is a large decrease in the third ligand 
binding energy. The trend for Fe(CO)„' is also different from 
ftat for the nfutral Fe(CO)„, where there is a dramatic increase 
in the second CO binding energy followed by a decrease for 
the third. 3 However for both Fe(CO)„- and Fe(CO)„ the third 
binding energy is similar to the fourth. In this work we consider 
the Fe(CO)„~ systems with n — 1—4. Our goals are o 
understand the nature of the bonding and to explain the trends 
in the experimental CO binding energies. 

In addition to the binding energies for all four CO m° lecu ^ s . 
there is detailed information about FeCO . Villalta and 
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II. Methods 

The Fe basis set is a [8s 4p 3d] contraction of the (14s 9p 
5d) primitive set developed by Wachters. 12 The s and p spaces 
are contracted using contraction number 3, while the id space 
contracted (311). To this basis set two diffuse p functions are 
addethtiiese are the functions optimized by Wachters multiplied 
by 1 5 A diffuse d function 13 and an f polarization function 
(a = 1 339) are added. To describe the negative ion, , a diffuse 
and potion (ct(s) = 0.013963 and a(p) = 0.02092) are 
added The C and O basis sets are [4s 3p] contractions of the 
(9s 5p) primitive set optimized by Huzinaga. 1 The ^s spacers 
contracted (5211). A d polarization function is added 
exponents are 0.75 for carbon and 0.85 for oxygen. Only the 
pure spherical harmonic components of the basis functions are 
Esed m all calculations. We perform one calibration calculation 
usSg a large basis set. The C and O basis sets are the 
augmented-correlation-consistent polanzed valence quadmpe 
7 eta basis sets 15 without the diffuse g function. The Fe basis 
Ss“c£«p lOd 6f 4gV[(6 + 1)S (5 + 1J ;4d 3f 2g] atomic 
natural orbital basis set^ with a diffuse s (0.012) and p (0.009) 
function added. 

In the density functional theory (DFT) calculations we use a 
modification 12 of the original Becke hybrid functional of the 
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StoSi « SUmmary ° f SuCCessive Bind “S Energies, 


Fe-CO 

FeCO-CO 

Fe(CO),—CO 

Fe(CO),—CO 


expt" 


33.1 ±3.5 

34.9 ± 3.5 
41.5 ± 3.5 

40.9 ± 2.5 


prese nt work 

B3LYP CCSDfiy 2,t* pop 
0.55 


19.7 

34.4 

45.3 

384 


23.1 (24.7) 


0.59 

0.82 

0.74 
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Sng L B3LYf° ,nPUted Geometrical Parameters 


to0k?,c erl, ?K Wang - and S<1UlreS ' The ^lues have been converted 
to compute the zero-point 

comfMrted using the large basis set. < The dissoc.ation energy for' fc-“ 
a?o , g S .v C e°„ mPUted '° *" F(3d ’ 4s2) ' The «> pfpS ls 

form 

( 1 - A )E x sl *“ + AE W + BE x Bc cke + CE c lyp + 

o - QE c vwn 

where E, sia ' er is the Slater exchange, £,«= is the H artree-Fock 

ofBecke’"' £ part of the exchan ge functional 
and Pai- 20 /vwn ' h corre,atlon functional of Lee, Yang, 
and w ’ -u ls the correlation functional of Vosko, Wilk 
d Nusair, and A, B. and C are the coefficients determined 

heat^o^fomat 8 f thr h ee ' parameter flt ^ the experimental 
formation for his original hybnd functional The 

modified functional is denoted B3LYP and is used to optimize 

^ t0 COmpUte the Pre< l uenc 'es. The computed 

to^iTdnima COnfirm that structures correspond 

enelvofFern™ 126 * egeomet ^ Md compute the binding 
energy of FeCO using the coupled-cluster singles and double! 

approach including a perturbationai estimate of the triples 
examnons 23 [denoted CCSD(T)]. The CCSD(T) approach 
based on a spur-unrestricted self-consistent-field wave function 
Chtly the valence electrons are correlated In addition we 

an^h ^ b ' ndlng energ y usin £ the CCSD(T) 

ligeb^r'S ?GSDm ' “1 *? ^ ^ B3LYP ge0meUy - These 

shell C^SdTv appr!2. '? ““■t* reSInC,ed 0p - 

exinff? C ,l! at j 0nS u Cre perf0rmed usin § Gaussian 92/DFT, 25 
except for the large basis set CCSD(T) calculations which were 

94 26 ^ v ' sua hzation system 

in. Results and Discussion 

d KCOh-. and F«CO»- we consider bod, 

doublet and quartet states. For FeCO‘ the quartet slate is 
significantly lower than the doublet state. For FeiCO i - the 

s^ et 2^ d teal/ l0U ! le h S | ateS ^ Cl ° Se 10 Cnergy Wilh ** doubl et 
state 2.5 kcal/mol below the quartet state. This energy 

of ~r t0 definitively determ, ne the ground state 
ot re(CO>2 To compute the binding energies, we use the 
doublet state as it is lower in energy. For Fe(CO) 3 ~ the doublet 
state is considerably lower than the quartet state. Because low- 

have T r by addit, ° nal ,igands - Fe(CoT- mu"t 

have a doublet ground state, and therefore, we study only this 

fnJ!* COmputed binding energies, geometries, and harmonic 

nZeofTheT ^ ‘“ 3 We f,rSt discuss the 

ture of the bonding. The ground state of FeCO - is 4 Z - with 
an Fe occupation of 3d < 3 2 3d.T 4 sda : s P CT l . This state is derived 


FeCO - 4 I - C_, 
Fe(COb - ! n u D_* 
FefCOb - 4 n„ 
Fe(CO), - 2 Ai' Dy, 

Fe(CO) 4 -‘ 2 /t,c 3 u 

axial 

equatorial 


r e (Fe— C) 


r e (C— 0)“ 


1.801 (1.860) 
1.815 
1.805 
1.797 

1.794 

1.799 


1.177(1.194) 

1.180 

1.176 

1.175 

1.165 

1.167 


a For comparison the bond length of free CO is 1 130 and mil 
at the B3LYP and CCSD(T) levels of theory, 

ccs£rt“ 8 ,S 99 90 and the Fe - C -° ^8‘es Ire 180 4 *£ 

in A. ( P me ' erS g,ven ln Parentheses. The bond lengths are 

B^VT E to cm*™ 00 * 0 Fre< ' uendes Computed Using 

CO 


C-0 

2211 





bend(.T) 


FeCO - 



208 

Fe-C 

437 

c-o 

1864 

bend(jr u ) 

wag(^ 

Fe— C(a u ) 
C— 0(a g ) 


Fe(COb - 

4 n u a 



73 

393 

452 

1919 

bendlTf,) 

Wag(7T u ) 

wag^) 

91 

432 

455 

wag(jr g ) 

Fe-C(a e ) 

C-0(<7„) 

324 

444 

1866 

bend(.T u ) 

wag(jr g 0 

wag(jr u ') 

C-0(a u ) 


Fe(CO)c 

2 n,° 



53 

385 

436 

1900 

bend(7r„') 

wag(zr u ) 

Fe-C(a„) 

60 

402 

466 

wag^) 
Fe-C(a g ) 
C— O(o 6 ) 

323 

433 

1860 

bend(e') 

wag(a 2 ') 

wag(a 2 /4 ) 

C-0(a,') 


Fe(CO), - 



56 

366 

558 

1978 

bend(a 2 ") 
Fe— C(a',) 
wag(e') 

71 

434 

583 

wag(e") 

Fe-C(e') 

C-O(e') 

285 

506 

1875 

bend(e) 
wag(e) 
Fe— C(e) 
wag(ai) 
C-O(ai) 


Fe(CO) 4 

- b 



62 

341 

472 

587 

1940 

bend(e) 

wag(a 2 ) 

Fe-C(ai) 

wag(e) 

C-0(a,) 

80 

365 

492 

591 

2026 

bend(ai) 

Fe-C(a,) 

wag(e) 

C-O(e) 

86 

421 

505 

1925 


The calculations are performed as n„ or n„„ hence the degenerate 

Z^!f rc . t Spl,t and we 6 , “‘ b 0 * components; the second component is 
denoted with a pnme. * The Fe-C sketches and CO wags at 505 587 
and 591 cm ' are somewhat mixed in character. 

from a 3d 7 4s 2 Fe - configuration, but a 3d population of 6 61 
indicates that a mixing of 3d‘ character alsoZcurs as a result 
the large Fe to CO 2ji* donation. The 3<Lt* occupancy 

““ ““ ? d “ C0 ^ d0 "“™ ™» *>„do„ S 

n a C O bond length considerably longer than that in free CO 
Md longer than that in FeCO and a C-O harmonic frequencv 

™i! S CI " sma,,cr than that m free CO and 167 cm - ' 
smaller than that in FeCO at the same level of theory. 2 * Thus 
the addition of an electron to Fe increases the 3d to CO 2 Jt* 
donation relative to FeCO, as expected. The (Fe 4s)-(CO 5a) 

~*^"*~*. but 11,6 repulsion Educed 
by bmh sda and sp hybndization. The extra a electron in 

FeCO relative to FeCO results in an Fe-C distance that is 
larger for FeCO - than for FeCO. 

2 ? u ^ 7 d 4 P U St3teS ° f Fe(CO)2 ~ have D»* symmetry and 
have an Fe 3d 4s‘ occupation. That is, the bonding is Fe°(CO)C 

^ e p thC CleCtr ° n ° n the (C0)2 subunit is in ^ 

Sh rh o Z 6 6 ° CCUpa,i0n is 3d<5 2 3d7T 4 sda 2 ,n both states. 

'* h ^ e'ectron high-spin-coupled to the open-shell 3d 
!** T0 "f ln quartet state and low-spin-coupled in the doublet 

to FeCn“ S b ? ndlng " lechanism is somewhat different from that 
FeCO In FeCO there are three Fe a valence electrons 
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and the Fe— CO repulsion is reduced by sdo and sp hybridiza- 
tion. Clearly, this mechanism cannot occur for Fe(CO)2 , if 
FelCOb is linear, sp hybridization is lost, while if it is bent, 
as has been found in other ML 2 cases, sda hybridization is lost. 
Thus, in Fe(CO) 2 ' one of these a electrons is transferred to the 
CO 2ji* orbital. We should note that this can also be viewed 
as a two-step process, where the Fe“ is promoted to a 3d 7 4s>- 
4p' occupation and the 4p electron is donated to the CO 2 ,t* 
orbitals. 

The ground state of Fe(CO) 3 " is 2 A', with D 3h symmetry. 
The bonding is perhaps easiest to view as being derived from 
a 3d 9 occupation of Fe“. Assume that the Fe atom and CO 
molecules are in the xy plane, with the Fe atom at the origin 
and one of the CO molecules along the y axis. To maximize 
the metal to CO 2 jt* donation, the 3d K and 3d y: orbitals are 
doubly occupied. The 3d„ and 3dj-y> orbitals hybridize with 
the Fe 4p, and 4p v orbitals, respectively, to maximize the metal 
to CO 2 n* donation with the CO along the y axis and to 
minimize the repulsion with two 5a orbitals of the other two 
CO molecules. This effect is clearly illustrated in a plot of the 
4pjr-3dJt hybrid orbital-see Figure 1. The orbital 

has the same symmetry as the symmetric combination of the 
CO 5a orbitals and is singly occupied to minimize the Fe— CO 
repulsion. We should note that because of the pd hybridization, 
the bonding in this molecule can also be viewed as being derived 
from the 3d 8 4p‘ or 3d 7 4p 2 occupations. As is clear from this 
analysis of the bonding, the Du structure is very favorable as 
it minimizes the ligand-ligand repulsion, and the polarization 
of the 3d orbitals that maximizes the donation to the CO 2n* 
orbital also reduces the Fe-CO repulsion. 

The lowest state of Fe(CO) 4 ' is 2 Ai and it has a C 3 » 
symmetry. This system is derived from FefCOb by adding 
an axial CO molecule along the z axis. The three equatorial 
CO molecules that are in the xy plane in Fe(CO);, bend out of 
the plane by 10°. The small bend retains the very favorable 
bonding described for Fe(CO) 3 '. The major difference in the 
bonding between Fe(CO) 3 ' and Fe(CO) 4 - is the singly occupied 
orbital. To reduce the repulsion with the axial CO, the 
3d2b-r--r hybridizes with the Fe 4p ; orbital. The unpaired 
electron is located opposite to the axial CO and between the 
three equatorial CO molecules. 

On the basis of their experiments, Villalta and Leopold 9 
predicted the equilibrium bond lengths of the X 3 2 and a 5 Z 
states of FeCO relative to those in FeCO'. In the analysis of 
the experimental data, one must make an assumption about 
relative change in the Fe-C and C-O bond lengths when the 
electron is detached. They performed the analysis making two 
assumptions, that both bond lengths could change in the same 
direction or in the opposite directions, but concluded that 
changes in the opposite direction were most likely. Our 
calculations show that the changes in the C-0 and Fe-C bond 
lengths relative to FeCO' do occur in opposite direction for 
the 5 2“ state; the C-0 bond length decreases and the Fe-C 
bond length increases. However, for the 3 2' state we find that 
both bond lengths contract. That is, the C-0 length decreases 
because there is less Fe to CO 2 jt* donation and the Fe-C 
distance decreases because the removal of the a electron 
decreases the Fe-CO repulsion. On the basis of our calcula- 
tions, we compare with experimental analysis that agree with 
the signs of our computed changes in bond length. The B3LYP 
Fe-C bond length of 1.801 A for FeCO' is between die B3LYP 
values of 1.770 and 1.912 A 28 for the 3 I and 5 2 states of 
FeCO, respectively, but closer to the 3 2 state. The experi- 
mental differences for the 3 X' and 5 Z' states are -0.05 ± 0.02 
and +0.07 ± 0.02 A, respectively, which are in reasonable 



Figure 1. Isosurface of the pd hybrid orbital at cut-off values of ±0.12 
au. The Fe(CO) 3 _ molecule is in the plane of the paper. 


agreement with the computed values. The experimental CO 
contraction, relative to FeCO", is 0.03 ± 0.01 A for both states 
which is in good agreement with the B3LYP contractions of 
0.028 and 0.023 A for the 3 2' and 5 2' states, respectively. The 
CCSD(T) results for FeCCT and for the and 5 Z states of 
FeCO 28 - 29 are consistent with those obtained using the B3LYP 
approach Castro, Salahub, and Fournier" also report similar 
trends for the Fe-C and the C-O bond lengths using the 
LCGTO-DF approach. The B3LYP Fe-C stretching and the 
bending harmonic frequencies are in good agreement with 
experiment. The C-O stretching frequency of FeCO has not 
been measured. Castro, Salahub, and Fournier predict a C-O 
frequency of 1831 cm' 1 which is close to our B3LYP value of 
1864 cm' 1 . Thus, excluding the binding energy, the results of 
Castro, Salahub, and Fournier agree with those obtained in this 
work. 


Theoretical Study of Fe(CO) n ' 

We now consider the trends in the CO binding energies. We 
first note that bonding in the positive ion is determined mostly 
by minimizing the Fe-CO repulsion, which is very different 
from the neuUal or negative ion, where Fe donation to the CO 
2ji* orbital is the most important factor in describing the 
bonding. Therefore, the difference in the trends in the CO 
binding energies for Fe(CO)„- and Fe(CO)„+ is not too 
surprising and is easy to explain based on this work and that 
presented previously. 2 The present calculations explain the 
difference between the negative ion and neutral, as well as the 
trend on the negative ions. To aid in the discussion we note 
that the CO binding energies (at 298 K) for the neutral systems 
are 42.9 ± 3.5, 36.7 ± 3.5, 29.1 ± 5.8. and 27.9 ± 8.8 kcal/ 
mol, where the first CO binding energy is given with respect 
to Fe F(3d 4s ). as this asymptote is to which it dissociates 
and is the most consistent with the dissociation for FeCO - . 
Based on the corrections that we computed for Fe(CO)„ - . these 
Fe(CO)„ values would probably be reduced by 0.4— 0.8 kcal / 
mol if corrected to 0 K. 

The computed Fe(CO) n - D 0 values are compared with 
experiment 3 in Table 1. The trend in the computed D 0 values 
agrees with experiment. Namely, the CO binding energies 
increase until the third, and then there is a small decrease for 
the fourth. This trend and the differences between the negative 
ion and neutral can be understood in terms of the bonding 
described above. For FeCO - , the extra electron results in more 
Fe to 2/T* donation, which enhances the bonding relative to 
that in FeCO; however, the extra electron increases the a 
repulsion. The sda and sp polarizations are efficient at reducing 
the repulsion for both the ion and neutral, but the extra a 
repulsion results in the negative ion being about 10 kcal/mol 
less strongly bound. In Fe(CO) 2 - the promotion of the electron 
from the 4s to the 4p orbital enhances the bonding, by reducing 
the a repulsion and increasing the donation to the 2 it* orbital 
but the binding energy is naturally reduced by the cost of this 
promotion. The similar first and second binding energies for 
Fe(CO)„ indicate that the promotion energy is very similar in 
magnitude to the enhanced binding. However, the binding 
energy of Fe(CO) 2 - is smaller than that in Fe(CO) 2 because 
Fe(CO)2 does not have to pay the s to p promotion energy. For 
Fe(CO) 3 the promotion energy is now shared by three ligands 
and the binding energy is larger than in Fe(CO) 2 -. The binding 
energy m Fe(CO) 3 is significantly larger than that in Fe(CO) 3 
because the extra electron in the negative ion allows for a much 
larger donation to the 2jz* orbitals. The binding energy of 
Fe(CO) 4 is slightly smaller than that of Fe(CO) 3 - . This is 
not too surprising given the very similar binding in these two 
systems. The extra CO results in a larger ligand-ligand 
repulsion and slightly smaller 2r* donation per CO This results 
in a decrease in the binding energy even though the promotion 
energy is now shared by four CO molecules. The extra electron 
again results in the negative ion being significantly more 
strongly bound than the neutral. 

In addition to the trend in the computed D 0 values being 
consistent with experiment, the computed B3LYP values for 
FefCOh , Fe(CO) 3 ~, and Fe(CO) 4 - agree well with experiment- 
the computed values for Fe(CO) 2 - and Fe(C0) 4 - fall within 
the experimental error bars and Fe(CO) 3 “ is only 0.3 kcal/mol 
larger than the experimental upper bound. Unlike these values 
the computed B3LYP result for FeCO" is significantly smaller 
““ * e ex Penmental result. The results obtained using the 
CCSD(T) approach are larger than the B3LYP result but are 
still much smaller than that from experiment. It is very difficult 
to assign all of this difference to an error in experiment because 
the CCSD(T) result 29 for FeCO is about 2 kcal/mol smaller than 
the most accurate experimental value, 9 which is derived from 
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i : Con, iwrison of the Computed Fe Electron 
Affinity (kcal/mol) with Experim ent 

B3LYP CCSD(T) expr 


Fe 4 F(3d 7 4s 2 )— Fe 5 D(3d 6 4s 2 ) 
error 

a Reference 30. 


14.5 

10.7 


- 7.8 

- 11.6 


3.5 ±0.1 


the experimental FeCO - Z) 298 value. Both the B3LYP and 
CCSD(T) agree 28 with experiment for FeCO + , while the B3LYP 
result 28 is too small for FeCO and the CCSD(T) agrees 29 with 
experiment, and both are too small for FeCO - . This suggests 
that as the Fe to 2jz* donation increases, the system becomes 
harder to describe using a single reference-based approach. We 
should also note that we have found that the error in the B3LYP 
approach is often the largest for the first binding energy and 
that part of this error arises from a description of the metal atom. 
In Table 4 we compare the computed and experimental 30 
electron affinity (EA) of Fe. The B3LYP value is significantly 
too large because the B3LYP is biased in favor of metal 
occupations with a larger number of 3d electrons. The Fe 
population in FeCO - is between Fe and Fe - . and the 3d 
populations is between six and seven. Thus, a second method 
of computing the binding energy would be to dissociate to CO 
+ Fe D(3d 6 4s 2 ) + e - and correct this to the CO + Fe - 4 F(3d 7 - 
4s ) asymptote using the experimental EA. This results in an 
FeCO B3LYP binding energy of 37.0 kcal/mol. This argument 
is too simplistic, however, as applying this approach would make 
the CCSD(T) result in worse agreement with experiment. Thus, 
while this error in the description of Fe - is probably the origin 
of part of the error in the B3LYP approach, it cannot be used 
to compute a more accurate value. Despite the small D 0 value 
for FeCO , the experiments of Villalta and Leopold 9 confirm 
that we have correctly described the nature of the bonding in 
FeCO and thus that our analyses of the bonding and trends 
are correct. 

IV. Conclusions 

The computed changes in geometry between FeCO - and the 
X X and a 5 X states of FeCO are in reasonable agreement with 
experiment, as are the computed FeCO~ frequencies. However, 
the computed Fe - -CO binding energy is too small. The 
second, third, and fourth CO binding energies are in good 
agreement with experiment. For FeCO - the calculations show 
that sda and sp hybridizations reduce the Fe-CO repulsion 
With the addition of the second CO molecule, there is a 
promotion of one of these Fe a valence electrons to the CO 
2rt* orbital. This reduces the o repulsion and increases the 7 i 
bonding. However, the first and second CO binding energies 
are very similar because of the cost of this promotion. In 
re(CO)3 and Fe(CO) 4 ~ the cost of this promotion is shared 
by more ligands, resulting in third and fourth CO binding 
energies that are significantly larger than the first two The 
bonding in Fe(CO) 3 ~ and Fe(CO) 4 - is very similar and the 
smaller 2 ji* donation per CO and the larger ligand— ligand 
repulsion results in a slightly smaller binding energy for the 
fourth CO than for the third. The calculations also explain the 
difference in the trends in the binding energies for the neutral 
and negative ion. 
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